Marginal scaling scenario and analytic results for a glassy compaction model 
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A diffusion-deposition model for glassy dynamics in compacting granular systems is treated by 
time scaling and by a method that provides the exact asymptotic (long time) behavior. The results 
include Vogel-Fulcher dependence of rates on density, inverse logarithmic time decay of densities, 
exponential distribution of decay times and broadening of noise spectrum. These are all in broad 
agreement with experiments. The main characteristics result from a marginal rescaling in time of 
the control parameter (density); this is argued to be generic for glassy systems. 
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Glassy dynamics occurs with similar characteristics in 
a remarkably diverse range of systems Q. This letter 
attempts to give reasons for this similarity. 
It begins by considering an idealized model of granular 
materials, which is perhaps the simplest class of real sys- 
tems showing glassy behavior. The model is treated by 
approximate scaling and then by asymptotically exact 
methods. The behavior is similar to that found in simula- 
tions and in the Chicago experiments [| ^, M, ^| on real 
granular systems. Because we expect diverging charac- 
teristic times (while no diverging length has yet been seen 
in glassy systems) the scaling procedure rescales time, t. 
The consequent scaling equation for the control parame- 
ter, in this case the density p, has a particular (marginal) 
form. We show at the end of the letter that this leads di- 
rectly to the best known characteristic features of glassy 
dynamics. We also argue at the end that this marginal 
form arises when time scales are excessively sensitive to 
changes of the control parameter, as is the case in glassy 
dynamics in general, so arriving at a universal scaling 
scenario. 

All the following characteristics of glassy behavior Q |t]] 
are shared by super cooled molecular glass-forming liq- 
uids, structural glasses, foams, colloids and even (shaken) 
granular systems: (i) extreme slowing of rates with typ- 
ically a Vogel-Fulcher dependence on control parameters 
(temperature, or density, etc); (ii) slow decay of correla- 
tion functions and of the control parameter itself, often 
fitted to stretched exponential or inverse logarithmic time 
dependence; {Hi) associated aging phenomena; (iv) char- 
acteristic broadening of noise spectra. 
This has raised questions of universality, and it suggests 
that, despite obvious differences in fundamental mecha- 
nisms, something generic might already be present and 
discernable in simple models. Granular materials show 
all the standard characteristics when compacting under 
shaking, even though thermal aspects play no role ||. 
The simplicity of models (e.g. exclusion models) appro- 
priate for shaken granular systems [§| |To[ [□], [l2|, [l3|, 
|l5f makes them ideal candidates for such investigation, 
and we consider one such model here. 
Previous investigations have shown g, [|, [l6) how a sim- 



ple one-dimensional deposition-evaporation model (con- 
tinuum car parking) can yield characteristic behavior. 
Even the noise spectrum (obtained through simulations 
on the model Q ) agrees well with experimental observa- 
tions §,§. 

The model considered here is closely related to the con- 
tinuum car parking model and consists of unit sized 
blocks performing a random walk on a ring of size L. 
They only interact through hard core potentials and as 
soon as a gap of unit size opens up between two adjacent 
blocks, an additional block will be deposited in the gap 
and the diffusion will continue (see Fig. [l]). 



i 

n 



in n ni i 1 in nn n n 



FIG. 1: Diffusion-deposition process on a ring 

This model allows rearrangement (diffusion) of the 
grains to enable another grain to move (deposit). The 
model mimics the actual mechanism of rattling within 
and escape from/entry into rearranging cages formed by 
surrounding grains. Under compaction, typical free vol- 
ume elements become smaller, drastically reducing the 
probability of finding one big enough to accommodate a 
block, and hence lengthening time scales. 

A qualitative scaling description of this effect is as fol- 
lows. Consider elementary walkers each carrying a gap of 
size e = 1 — p, which is the average free volume per block. 
For the formation of a gap of unit size, n = 1/e such ele- 
mentary walkers have to coincide. A typical initial state 
has on average one elementary walker on each lattice site. 
Since walkers can move a distance of order \fi in time t, 
by t ~ n 2 walkers from n sites around a particular site 
can have accumulated there, and do so with approximate 
probability per unit time 



(t-V2 CX p(-n 2 /t)) n 



(1) 



It is straight forward to show that this result applies 
in any dimension. This rate of unit gap formation is 
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equivalent to a characteristic inverse time t _1 for depo- 
sition/motion from or into a cage. Equation (|l|) implies 
that under rescaling of r by a dilatation factor b, n goes 
to n' where n' n = bn n . This equality can be written 
in terms of the free volume density parameter e as the 
scaling equation 

e' = e + Ae 2 \nb+---. (2) 

In this treatment A ~ 1/lne. However, the exact discus- 
sion of the asymptotic behavior obtained below, implies 
that the proper form of (|^) has A a constant. In either 
form, the scaling equation has a marginal form, in which 
to leading order the parameter does not rescale. This 
marginality can also be shown to occur in the continuous 
car parking model and we argue below that this is a key 
feature of all adequate models of glassy behavior, having 
its physical origin in the exceptional slowing in such sys- 
tems. 

First we turn to the full analytic treatment of the model. 
At some instant there are n blocks on the ring of size L. 
We choose to consider the gaps between adjacent blocks, 
Aj, i = 1, . . . , n, as the dynamical variables. The vector 
A n = (Ai, . . . , A„) then performs a random walk on the 
hyper-surface (See Fig. ||) 




FIG. 2: A 3 (t) tracing out a random walk trajectory on its, 
being reflected once at 7r 3 |A 3 =o and then passing through the 
portal at 7T3|ai=i 

tion process we impose the boundary conditions: n n \ Ai=o 
reflecting (corresponding to two blocks bouncing off each 
other) and 7r„|A i= i a portal that will transfer the system 
to ir n+ i |A i= A i+ i=o (corresponding to a gap of unit size 
forming and a block being deposited). 
In what follows, unless stated otherwise, ~ means 
"asymptotically correct in the large system, high density 
limit" , and the quantity p (or equivalently e = 1 — p) will 
be used as our control parameter. The time evolution of 



e(t) is governed by 

m--Tkr } => (m)=-L- 1 (l/r(t)). (3) 

Here r(t) is the time between deposition events at time t 
(i.e. the total time spent on Tr n u\), for the specific realiza- 
tion considered. In order to estimate (l/r(t)) we define 
the cross-section <r(e) as the ratio between the portal area 
and the total boundary area (each (n — 2)-dimensional) 
on each hyper-surface 7r„. This will allow us to find the 
probability of hitting a portal, given that we hit a bound- 
ary. Through direct volume consideration one gets, after 
some lengthy analysis, the asymptotically exact result 

(7(e)~exp(l-l/e). (4) 

The same analysis tells us that the distances between 
boundaries of ttl p scale as e. Therefore if we let Af(e) 
be the average time it takes for A Lp to go between 
boundaries, we have Af(e) ~ fee 2 , where A: is a con- 
stant inversely proportional to the diffusion constant of 
the blocks. Define, with a slight abuse of notation, r(e) to 
be the time between deposition events at a given density. 
Since the portal cross section is very small for high densi- 
ties we expect the random walk to be essentially ergodic 
between the deposition events. Therefore we can get the 
distribution for r(e) through considering repeated bounc- 
ing off the reflecting boundaries (blocks colliding) until 
finally hitting a portal (deposition event). The proba- 
bility of hitting a portal on the k:th contact with the 
boundary is then (1 — cr) fc_1 (T, so the probability distri- 
bution for r(e) is 

P(T{e)=x)~{l-v(e))*W>- 1 v(e), (5) 

where a; is a multiple of Af(e). In the limit e -> (||) 
can be considered to be a continuous exponential distri- 
bution. By increasing L we can, for any fixed density, 
have an arbitrary number of essentially independent de- 
position events in any finite time interval. From the law 
of large numbers we thus know that the e(i)-distribution 
can be made arbitrarily sharp in the large system limit. 
Therefore r(t) ~ r((e(f))). To solve |§ for (e(t)) we 
need to calculate the average deposition rate at time t, 
(1/t(*)) - (l/r((e(t)))). This is easily done using (|), 

(W))~^ F exp(-l/(e(t))), (6) 
which inserted into (j|) gives the asymptotic solution 

~ In(etfkL) ~ Int CD 

for large t. The asymptotic solution (]?]) is not valid at 
densities low enough to be suitable for comparison with 
practical simulations. Since the cross section (Q) becomes 
small (i.e. ergodicity holds between deposition events) 
already at quite moderate densities we instead compare 
a numerical integration of (pi) to simulation data (see 
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FIG. 3: Least square fit, with respect to the constant fc in (|6j) , 
of the numerical integration of (pi), compared to data from a 
simulation of the evolving density using parallel updates on a 
ring of size L = 10, 000. 



Fig. |). 

To examine the power spectrum we need to include the 
effects of fluctuations. Since ergodicity holds on tt, in the 
densely packed limit, we can write the Master equation 
for the probability distribution of p at time t, P(p, t), as, 



1/L,t) 



r( P ) 



P( P ,t). 



This is easily solved by implementing Laplace transform 
techniques on the time variable. For the initial condition 
p(0) = pi and in the continuum limit the solution is (to 
the highest order of L in the exponential) 



P(p,t\pi) ~ r(p)J 2vT 2 ( P , Pl ) ex P 



where p is the continuum probability density and 



(Sp(t A +t)Sp(t A )) 



2L l+t/t A ' 



(t > 0). 



(9) 



In the last step we have used our particular form of r(p) 
and assumed that we are in a region where e ~ e A ~ £/, 
but where t(c),t(c a ) ^> r(e/). The above result shows 
that our model exhibits "simple aging" in the sense that 
the two-time correlation function is only a function of 
tit a- The density-density response coefficient, corre- 
sponding to perturbing the density at time t A and then 
measuring the response in the density at time t A + t, is 
proportional to (^|) (ignoring an additive constant) (||). 
This gives the simplest form of the fluctuation dissipa- 
tion violation investigated in jl7], |l8| . We now define the 



(complex) power-spectrum for the aged system as 

at a 

S a (uj 7 t A ) = J dtcxp(iuit) (Sp(t A + t)Sp(t A )) . 
o 

Letting E\ (z) denote the exponential integral, this can 
be written in the scaling form 

S a (uj,t A ) =t A ^s a (t A uj), 

s a (x) = e lx (E 1 (ix)-E 1 (ix{a + l))). 

Since the rates in our system are of the order of the in- 
verse aging time (i.e. r(e A ) ~ t A ), a fairly large a must 
be considered in order to sample the effects of aging. 
Therefore there are three asymptotic regimes (see Fig. ||) ; 



A) l/a < \x\ < 1, s a (x) 

B) 1< \x\ < a, s a (x) 

C) a <C s a (x) 



ln(a + 1) — (a — hia)ix 
l i 



sin(ax) 
x(a + l) 



(10) 



The reason for the lower cut-off in x is that the values 
of to < l/at A are unphysical. Further we know that on 
short time-scales the system is time translation invariant 
and therefor the real part of S a corresponds to the "or- 
dinary" power-spectra in this limit. We naively stretch 
this analogy to the whole range of possible uj (and note 
that the the imaginary part of S a stays close to the real 
part of S a for the longer time scales (see Fig. ^)). 



T n (p,Pi) — J dp'T n (p'). In order to keep the density $ a 
pi 

evolution rate finite in the continuum limit r has been 
rescaled as r — + Lt. Since r(p) diverges strongly as 
p — * 1, the approximations used to arrive at (ph break 
down, at any finite L, in the dense limit e < 1/vL- How- 
ever, for large L the average of p is independent of L, 
while the variance can be made as small as we please. So 
for large L, matching of (||) to a Gaussian in p gives the 
correct distribution for large L. Using the two variable 
Gaussian approximation the two-time density correlation 
function can be calculated for a system which starts at 
density pi and is aged for a time t A from whereon mea- 
surements are performed: 
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FIG. 4: Log-log-plot of the real and imaginary parts of the 
power spectrum scaling function s a (x) versus x = ujtA for a = 
"measuring time" /t a ~ 200. The as ymp totic regions and the 
asymptotic behavior (see equations (llOj)) are indicated. 

Since our model only allows deposition of blocks 
(grains) it corresponds to the limit of very weak tapping 
in ||. Taking this into account a qualitative agreement 
can be seen between our analytical results (see Fig. ^) 
and the experimental results presented in Further, 
assuming that the density increases with the depth we 
can interpret this as an increase of the effective aging 
time as measurements are done closer to the bottom of 
the container. Doing this we see that the bulk of the 
power spectra is moved towards lower frequencies while 
the maximum height seems to increase with aging. This 
is in full accord with the scaling form of S a . 

We now discuss further the nature of time scaling in 
such models and the question of universality. Under 
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rescaling of time (t — > t' — bt), the general renormal- 
ization group transformation [ pl| of a control parameter 
p, such as the density (or temperature) can be written 
as p — > p' = Rb(p)- If there is a divergence of the char- 
acteristic time scale r at p = p* , then p* is a fixed point 
of the transformation. And linearizing about the fixed 
point using e = p — p* gives in general 

e' = X b e + 0(e 2 ) (11) 

where (from the semi-group property RbRb> = Rw) the 
eigenvalue Xb (if non-zero) has to have a standard power- 
law dependence on the time dilation factor b, of the form 
A 6 = b 1 ^, with y > 0. This leads to r - \e\~ y , which is 
ordinary critical slowing. But the essential characteristic 
of glassy behavior in general is that the slowing of char- 
acteristic rates with the control parameter is excessive, 
more extreme than with ordinary critical slowing. Then 
the control parameter has to change only marginally to 
stretch the time scales significantly. That is the case 
when, in (pr|), Xb = 1. Then higher order terms in the 
expansion in e have to be retained, and ( |ll|) becomes of 
the marginal form 

e' = e + ae 1+Q ln6+ . . . , a>l 

where the 6-dependence is again set by the semi-group 
property pl| . As was seen earlier, this form, with a = 1, 
applies to the diffusion-deposition model (and also the 
continuum car parking model, etc). Eliminating b be- 
tween the scaling equations for e and for time gives a 
generalized Vogel-Fulcher form: 

T ~ ex P (ssf) ■ 

This specific equation is a direct result of the qualitative 
"stronger than ordinary critical slowing" essential 
characteristic of glassy dynamics. It suggests that the 
common Vogel-Fulcher form (a = 1) is one of various 
generic universal forms distinguished by the value of 
a. The second commonest form would then be the 
exponential-inverse-squared form (a — 2), which has 
appeared in a glassy model with constrained dynam- 
icsf||J. 

A second characteristic of glassy systems, including 
the models discussed here, and responsible for meta 
stability, is that the control parameter is one whose mean 
and local fluctuations evolve with the slowed internal 
dynamics. Its evolution towards the fixed point value 
then becomes of form oc (lnf) -1 /" when the generalized 
Vogel-Fulcher form applies. This scenario applies, with 
a = 1 , to the analytic treatment of the density evolution 
of the diffusion-deposition model, equation (j^). The 
slow evolution of the control parameter is in general 
related to the broadening log frequency scale of the noise 
spectrum. This time-scaling scenario seems equally 
applicable for 'thermal' glasses. It does not require a 



diverging length in glasses, nor a real static equilibrium 
transition. It does not so far explain possible/universal 
power laws in the noise spectrum, nor the non universal 
differences which distinguishes various types of "glass" . 
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